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In PrOs4Sbi2, the lowest-lying singlet and triplet states in a Pr 4/^ configuration hybridize 
with conduction electrons having local a„ and tu point-group symmetries. It is shown that for 
an attractive triplet pairing interaction, the orbital degrees of freedom of the tu component are 
important. In addition, the Th point-group symmetry characteristic of skutterudites plays an 
important role in stabilizing triplet superconductivity. 
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The skutterudite PrOs4Sbi2 is a recently found new 
heavy fermion superconductor.^ It is reported that the 
time reversal symmetry is broken in the superconduct- 
ing state, ^ and that there are multiple phases depending 
on a magnetic field. These results together with the ab- 
sence of the Hebel-Slichter peak in nuclear quadrupole 
resonance (NQR) measurements^ indicate that uncon- 
ventional superconductivity is realized in PrOs4Sbi2. In 
theoretical studies, an attempt to confirm these experi- 
mental results of unconventional superconductivity has 
been carried out.^"* 

The Pr'^+ ion in PrOs4Sbi2 has a 4/^ configuration in a 
Th (without fourfold axes in a cubic crystal) point-group 
crystal field.^ The Fi singlet ground state is realized 
with low-lying r4 ' triplet excited states. This crystal- 
field configuration can explain naturally the magnetic 
field-induced antiferro (AF)-type order above 4.5 T,^^"-'^'^ 
where the superconductivity disappears. This is the most 
important feature compared to the reference supercon- 
ductors LaOs4Sbi2 and PrRu4Sbi2.""^^ Thus, we be- 
lieve that the existence of the low-lying triplet state is 
important for superconductivity. 

A recent neutron scattering experiment revealed that 
the low-lying excitation mode has a weak dispersion re- 
lation (exciton). In our previous study, we investi- 
gated exciton-mediated superconductivity by taking only 
an effective magnetic exchange interaction between the 
Pr triplet states and conduction electrons. In this case, 
we found an attractive interaction for d-wave supercon- 
ductivity. 

The observed exciton dispersion indicates that there 
are intersite multipolc interactions between the Pr ions. 
In the field-induced ordered phase, the following exper- 
imental results are reported: (1) Only a small dipole 
moment is observed. (2) The critical field exhibits a 
marked anisotropy depending on the field direction. 
(3) An inelastic neutron scattering experiment revealed 
that the intensity of the exciton decreases towards the 
zone boundary. The above results support a nonmag- 
netic multipole (quadrupole) interaction. ^^"^^ To derive 
such a substantial multipole interaction via a conduc- 
tion electron system, we have to take the orbital degrees 
of freedom into account in the exchange interaction be- 
tween the Pr 4/^ and conduction electrons.^'' Their cou- 



pling should contribute to the mass enhancement of the 
conduction electrons.^' 

By neutron scattering measurements, evidence of an 
AF multipole intersite interaction was also found. It 
is known that AF magnetic fiuctuations tend to stabi- 
lize singlet pairing superconducting states. If triplet su- 
perconductivity is stabilized, orbital degrees of freedom 
are important, as discussed in the case of the p-wave 
superconductivity of Sr2Ru04.^^ We apply this idea to 
/-electron systems in a different manner and pursue a 
stability condition for the triplet superconductivity. To 
determine the role of orbitals, PrOs4Sbi2 is a good ex- 
ample, since its orbital is well defined by the localized 
nature of the /-electron. In fact, recent NMR and /iSR 
measurements reported evidence of odd parity supercon- 
ductivity.^'''^* 

As mentioned above, the ground state in Pr 4/^ is |Fi) 
accompanied by low- lying triplet excitations |F4^'^) (n = 
1,2,3). In the Th symmetry, the |F4„) states are linear 
combinations of F5 and F4 wavefunctions in Oh-^''^^ 

\r^2) = VT^K)+d\n) (1) 

Here, n=l,2,3, and d is the Th crystal-field parameter de- 
scribing the deviation from the Oh symmetry (d — 01 
\d\ = 1). The Pr ion is at the center in the cage of Sbi2. 
The 4/2 state hybrid izes with conduction electrons hav- 
ing local a„ and i„ symmetries, which can be formed 
from the molecular orbitals of Sbi2.2^ In this paper, we 
treat a single dominant conduction band with strong 
hybridization with the Pr /^ state. The a„ electrons hy- 
bridize with the / Fy states of Pr directly and also with 
the / Fg states via nonlocal mixing with the <„ molecular 
orbital of Sbi2. 

Let us begin with the following nonlocal exchange in- 
teraction between the 4/^ state of the zth Pr and the 
conduction electrons:^'' 

^e.= ^^^E^'' (2) 
k k' 

X Y.iS^^4ik, k') + Q.A9^:(fe, k')] ^L'^.*^ fcV'- 

A/i (T(t' 

Here, J = VT + 20fFjo with Jo as the coupling con- 
stant for d — Q. Six and Qix are local magnetic and 
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nonmagnetic operators, respectively, which couple IFi) 
with |r4^^). The A = +, — components are defined by 

^. = (|rf2^)(ri|) + (|ri)(|rif|), 

Q, = {\Tf^){T,\)-{\T,){Tf^\), 

s. = si = -y2(|ri)(rfi)|) + y2(|rf3^)(ri|), 
Q_ = Ql = -V2(|ri)(rif I) - V2(|rf3))(ri|). (3) 

In cq. (2), n takes /j, — 0, z, +, —, where cr" is a unit ma- 
trix and cr^ = (cr^±i(T^)/2. The momentum dependences 
of tlic /i components of the magnetic and nonmagnetic 
channels are described by s'^{k, k') and k'), respec- 
tively. The former is defined by 



stants for the magnetic and nonmagnetic interactions, re- 
spectively. From eq. (6), we derived the following Hamil- 
tonian for the crystal- field excitons:^^'^^ 



s'i{k,k') = i^A^iA 
.5 A, 



Ay^), s^(fe,fe') = 0. 



a+(fe, fe') = i--[—{A\, + AVi) + /3(A + A 
s~{k,k') = [s+(fc',fe)]*. 



a ■ 



(4) 



(5) 



VI + 20^2 VI + 20(i2 

We introduced the parameter A representing the am- 
plitude of the nonlocal mixing between the a„ (xyz) 
and tu (x,y,z) components of conduction electrons. The 
symmetry of the mixing is expressed by the func- 



tions A and A fg2' ^ ki = P 



Fkxy, where F 



kyz 



iF, 



Ak2 ^ 



sin(fca;/2) sin(fcy/2) cos{kz/2), 



kyz 



cos(fca;/2)sin(fcy/2)sin(fc2/2), F 



sm{kx/2) cos{ky/2) sm{kz/2). They reflect the symme- 
try of the bcc lattice. s'^{k,k') is obtained simply by 
s'^{k,k') = [sl^(fc', fc)]*. Here, /x = 0, and 
—fjb = 0,2;,—,-!-. For nonmagnetic channels, q'^{k,k') is 
obtained from s^(fc,fc') by replacing ~^ 
(-^ fc'i>-^ fc'2) for ^ = z,-, and (Afci,Afc2) ^ 
{-A fci,-A fc2) for A = +. 

Within the above exchange interaction (2), we ob- 
tained the following simple form for a multipole- 
multipole interaction for the nearest-neighbor Pr pairs:^^ 



H, = Y,{DsSi-Sj+D,Qi-Q^), 

{ij) 



(6) 



since x, y and z are identical in the Th symmetry. Here, 
Qi ■ Qj = -QizQjz + {Qi-Qj+ + Qi+Qj-)/2 due to the 
anti-Hermite nature of Qz defined by eq. (3). de- 
notes the summation over the nearest-neighbor Pr sites 
on the bcc lattice. and Dq in eq. (6) are coupling con- 



if 



k 



\yb ky + ^ kz^ kz), 



(7) 



3 A cx 

8'_{k, k') = i--[-^(A*fc, + A\,^) + [3{A fci + A fc-i)], 
st{k,k') = i^Ap{A fe2 + ^fe'2). 

sZ{k,k')^i^A^^{A k2+Au'2)- 

Here, a and (i are functions of the parameter d: 
VI - d? „ V21d 



where b is the bosonic operator for the a = x,y,z 
component exciton. ^ is the dispersion for exciton 
Ek ^ v/[A + {Ds + Dg)A fc]2 - [{D, - Dg)A k?, with 
A = 8 cos(fc2;/2) cos(fcj,/2) cos(fc2/2) for the bcc lattice. 
A (A ^ Ds,Dq) is the crystal field splitting from |Fi) 

to |rf ). 

For Pr-based systems, it was pointed out that the exci- 
ton plays an important role in mass enhancement'^" and 
superconductivity.'^^ In this paper, we derive the effective 
interaction iJeff between conduction electrons mediated 
by the exciton within the following second-order pertur- 
bation: Hqs = iJcx[l/(-Eo — i?o)]-ffex. Here, Hq is the 
unperturbed Hamiltonian for the Pr 4/2 and conduction 
electron systems. We can eliminate bosonic operators by 
taking expectation values at low temperatures. Since we 
are interested mainly in superconductivity, we only con- 
sider the fc'C^fc'C fcC_ fc type of interaction. We define 
the singlet and triplet pair operators by 



= -^^i(^la'(^ kaC- 



ka' 



(8) 



where, a = x,y, z. 

For singlet pairings, we obtain the effective Hamilto- 
nian 



h: 



kk' 



Us: 



A2 



^ i-P kxyP k'xy "I" ^ kyz^ k'yz ^ kzx^ k'zx)^\^ k' 



(9) 



where Usi = 10/9 and Us2 = 2{Ds + Dg/9). For 
a2/3 _ ^2 Q |-|^| ^ 1/8), dxy, dyz and dza;-waves are 
stabilized in the lowest order of 1/A. For a larger \d\, 
however, the lowest order term becomes repulsive. For 
|d| > 1/8, we have to consider the second-order term pro- 
portional to A y /A^. The even parity component of 
A also leads to d-wave states. In this case, the pair 
wavefunction, for instance, the d^jz-wave, has a form of 
[sin(fca;/2) sin(A;y/2) cos(fcz/2)][cos(A;x/2) cos(fcj,/2) cos(fcz/2)]. 
For a small Fermi wave number fcp, the d-wave is favor- 
able due to large cosine function values. However, we 
note that the coupling for the singlet pairing is strongly 
suppressed around \d\ ~ 1/8 (q:2/3 = 0^). 

For triplet pairings, an effective interaction similar 
to eq. (9) is derived. There is no first-order term in 
1/A for triplet pairings due to the even parity nature 
of the fmictions [F ^xyiF kyzi^ kzx)- The interaction is 
of the order of 1/A2, and proportional to products of 
{P kxyj P kyz ; F kzx) functions and the exciton dispersion 
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A . There is, for instance, the following term: 

\^ kyz^ k'yz^ k k kzx^ k zx'' fe'' fe' 

_l_ p p , +^tiz \ y ^ fc- fc' 
"t" feaj;-^ k'xy^ k k' ' ^2 

Since A y contains p-wave-typc terms such as 
sin(fca;/2) cos(fcy/2) cos(fcz/2), /-wave pairing functions 
are obtained by the product of this p- wave-type function 
and the d-wave-type functions (F ^.^y, F ^.y^, F j^^^). For 
a larger fcp, the /-wave pairing states take advantage 
against the d-wave, since they contain three sine func- 
tions with larger values. 

Next, we focus on such triplet pairing states. We ex- 
press the effective interaction for triplet pairings in Oh 
bases as follows: 

H^^'" = ^j^T.T.Urr4mr'{k'). (10) 

fcfc' rr' 

Here, V^s ^ (l2JA/7f(D, + D^/Q) is an effective 
coupling for triplet pairings. Both magnetic {Dg) and 
nonmagnetic multipole-multipole interactions con- 
tribute to the triplet superconductivity. 

5r(fc) = V^{k)i\ + V^{k)i\ + V^{k)i\ (11) 

is an operator for a triplet pair with the Oh represen- 
tation r. and Vp*(fc) represent the spin and orbital 
components, respectively. Urv represents an interaction 
connecting V and V . We introduce a d- vector by a lin- 
ear combination of d{k) = ^pArVr(fc)- Here, Ap is 
the amplitude of the order parameter for F. The order 
parameter Ar can be obtained by solving the following 
gap equation: 

Ar = -27rTA/-o V C/rr'(Vr'(fe)-Vr"(fe))fcAr". 

(12) 

Here, {■ ■ ■) k represents an integral over the Fermi sur- 
face. LUrn is a Matsubara frequency for fcrmions. A'o is the 
density of states at the Fermi energy. In the Th system, 
the gap equations can be reduced to the ri®r2, r3i®r32 
and F4 ® Fs types. We find that the transition temper- 
atures for these pairing states are similar. We focus on 
one of the pairing types (F4®F5) in this paper, since this 
three-dimensional representation is most appropriate for 
the observed experimental results. The V[?(fc) functions 
in cq. (11) for the F4 © F5 type are listed in Table I. The 
interaction Urr' in eq. (12) is given in Table II. The off- 
diagonal matrix elements are proportional to ap. They 
vanish in the Oh system. The eigenvalues of the matrix in 
eq. (12) determine the effective interaction corresponding 
to the obtained eigenmode. A negative (positive) eigen- 
value means an attractive (repulsive) interaction. The 
lowest eigenvalue leads to the highest transition temper- 
ature. The eigenvalue depends on the Th parameter d 
strongly via the off-diagonal components. 

Let us discuss effects of the Th symmetry on supercon- 
ductivity. As shown in Table II, the off-diagonal terms 
connect the lowest eigenstate with higher-lying eigen- 
value states. This effect is similar to the level repulsion 
between the states. As a result, the lowest eigenvalue is 



Table I. Orbital component of the triplet pair operator. Here, a:, 3/ 
and z represent x = shi{kx/2), y = s\i\(ky/2) and z = sin(fcz/2), 
respectively. In the same manner, x' = cos(fci/2), y' = cos{ky/2) 
and z' = cos(kz/2). Since the r4 © Fs type has a three- 
dimensional representation, there are other two components ob- 
tained by replacing (x-, y, z) (y, z, x) cyclicly. There are two 
basis functions (denoted by A and B) for the r4 type and three 
basis functions (denoted by A, B and C) for the Fs type. 
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Table II. Matrix elements of Urr'- m = (3/2)0^ - (l/2)/32, 
02 = -(l/6)a2 + (l/2)/32, as = (2/1/3)0/?, 04 = 3/2 and 
05 = (4/3) a^. 
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reduced compared with that in the absence of the off- 
diagonal terms, and the effective attractive interaction is 
enhanced as in Fig. 1. Therefore, the Th symmetry sta- 
bilizes the triplet superconductivity. In contrast to the 
triplet pairing states, there is no such effect for the d- 
wave states [see eq. (9)]. 

Next, we apply our result to explain the recent ex- 
perimental data. Zero-field /zSR measurements revealed 
an internal magnetic field emerging at the supercon- 
ducting transition temperature.^ This indicates that su- 
perconductivity with a broken time reversal symmetry 
is realized in PrOs4Sbi2, implying the multicomponent 
rsi ® F32 or F4 ® F5 type. For the F4 ® F5 type, there 
are threefold degenerate solutions. There are two types 
of combination with a broken time reversal symmetry. 
One is the dr^,^ (k) ± idp^^ (fc) type, and the other is the 
dr^^k) + udr^^ik) + u'^dr^^ik) type with u = e^''^''/^. 
These states with a broken time reversal symmetry are 
nonunitary states with a finite q{k) = id(fe)* x d(fc) vec- 
tor. 

Let us discuss the Th nature of the gap function 
^|d(fc)p — \q{k)\ for the nonunitary states. In Fig. 2(a), 
we show the gap function, for instance, d(fc) — d^^^ (fe) 
for a Th parameter a = (3 = l/V^ {d = 1/V22). Al- 
though the shape looks like a clover with four leaves, 
there are six point nodes along all the basal axes x, y 
and z. The fourfold symmetry around the z axis is bro- 
ken under the Th symmetry. The angle of one leaflike 
shape tends to be small towards the y direction, and the 



4 J. Phys. Soc. Jpn. 



Letter 



M. Matsumoto and M. Koga 



0.01 



-0.01 



-0.02 



-0.03 



Fig. 1. Tfi parameter d dependence of the effective interaction for 
r4 ffi Fs type in arbitrary units, fcp = 0.757r is considered here. 
Red and black hnes are the eigenvalues with and without the 
off-diagonal interaction terms, respectively. 





orbital fluctuations are expected. Thus, the parameter A 
can have a sufficiently large value of A > y/D/ A. In this 
case, the superconducting state discussed in this paper 
can be stabilized. The triplet pairing is promising when 
the singlet pairing interaction is suppressed [\d\ ^ 1/8 
in eq. (9)]. Since the spin-orbit interaction works in the 
exchange process via the 4/ orbitals, the spin of the con- 
duction electron is not conserved in the scattering with 
orbital exchange. This can give rise to attractive triplet 
channels even when the exciton has AF fluctuations. We 
also found that the Th symmetry stabilizes the triplet 
superconductivity and reflects the symmetry of the gap 
function. It is remarkable that the atomic nature of Pr 
is revealed by the superconductivity. 
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Fig. 2. Gap functions at = (3"^ = 1/2 (d = 1/V22) for fcp = 
O.TStt: (a) d(fc) = dr^^ik) and (b) d(fc) = dr^Jfc) ± idr^^{k). 



gap decreases strongly in the x direction. The other two 
states [dr^^{k) and dr^^ik)] can appear with a broken 
time reversal symmetry, since the three states are de- 
generate. As shown in Fig. 2(b), the dr^^ component is 
combined with dr^^ , and the gap function takes a value 
close to the y direction, which is larger than that for the 
single component of dr^^ ■ 

Multiple superconducting phases are reported by 
thermal conductivity measurements.'^ In the lower-field 
phase, it was reported that the gap function has a 
twofold symmetry around the z axis,^''^^ which may cor- 
respond to the nonunitary triplet pairing states d{k) — 
dr^^ (k) ± idr^^ (fe) discussed above. 

In our previous study, we found the order of JqD^/A^ 
for the coupling of the d-wave pairing interaction, where 
only an effective magnetic exchange interaction between 
the triplet excitations and conduction electrons was con- 
sidered.^^ Here, D (^ A) represents the width of the ex- 
citon dispersion. In contrast to this, we considered both 
the a„ and i„ local symmetries in the single conduction 
band, where the latter has orbital degrees of freedom. 
This model leads to attractive channels for triplet pair- 
ing with the order of A'^JqD/A'^. 

If the field-induced order is a quadrupolar type, strong 
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